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Abstract 

Propulsion by growing actin networks is a universal mechanism used in many different biological 
systems, ranging from the sheet-like lamellipodium of crawling animal cells to the actin comet tails 
induced by certain bacteria and viruses in order to move within their host cells. Although the core 
molecular machinery for actin network growth is well preserved in all of these cases, the geometry 
of the propelled obstacle varies considerably. During recent years, filament orientation distribution 
has emerged as an important observable characterizing the structure and dynamical state of the 
growing network. Here we derive several continuum equations for the orientation distribution of 
filaments growing behind stiff obstacles of various shapes and validate the predicted steady state 
orientation patterns by stochastic computer simulations based on discrete filaments. We use an 
ordinary differential equation approach to demonstrate that for flat obstacles of finite size, two 
fundamentally different orientation patterns peaked at either ±35 or +70/0/— 70 degrees exhibit 
mutually exclusive stability, in agreement with earlier results for flat obstacles of very large lateral 
extension. We calculate and validate phase diagrams as a function of model parameters and show 
how this approach can be extended to obstacles with piecewise straight contours. For curved 
obstacles, we arrive at a partial differential equation in the continuum limit, which again is in good 
agreement with the computer simulations. In all cases, we can identify the same two fundamentally 
different orientation patterns, but only within an appropriate reference frame, which is adjusted to 
the local orientation of the obstacle contour. Our results suggest that two fundamentally different 
network architectures compete with each other in growing actin networks, irrespective of obstacle 
geometry, and clarify how simulated and electron tomography data have to be analyzed for non-flat 
obstacle geometries. 
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I. INTRODUCTION 

The growth of actin networks is a generic propulsion mechanism occurring in a large 
variety of biological systems, ranging from the protruding lamellipodia of animal cells to the 
actin comet tails recruited by pathogens like the bacterium Listeria monocytogenes or the 
virus Vaccinia within their host cells. Due to its central importance, the molecular basis of 
this process is well preserved over a wide range of different species [Ij. Although a large num- 
ber of accessory proteins is known to be involved in actin dynamics on the cellular scale, in 
regions close to the leading edge the dynamics of network growth are determined by a small 
number of key reactions. Here the interplay between three fundamental processes deter- 
mines the structure of the growing network: polymerization of filamentous actin, branching 
mediated by the protein complex Arp2/3 and binding of capping proteins to the filament 
tip preventing further growth [21 |3]. The fact that these processes are highly conserved 
is impressively demonstrated by the observation that many intracellular pathogens rely on 
them for efficient infection and spread in the cytoplasm of their hosts [U-E] . Even function- 
alized plastic beads, rods and discs as well as lipid vesicles and oil droplet in purified protein 
solutions containing this minimal set of molecules have been shown to be propelled in vitro 
by polymerizing actin networks [7I4T2]. 

While the underlying molecular basis is very similar in all of these cases, the geometry 
of the different propelled objects is very different. For instance, the highly curved shape of 
relatively small, almost ellipsoidal pathogens like Listeria monocytogenes or Vaccinia virus 
is very different from the relatively flat shape of the leading edge of a crawling cells. The 
exact shape of the membrane in migrating cells is very dynamic, but certainly corrugated 
on different scales, thus in this case a flat obstacle shape can only be a flrst approximation 
and corrugated obstacle contours are also of large interest. In a growing actin network, new 
fllaments nucleate by branching off from existing mother fllaments at a characteristic angle 
around 70° set by the molecular geometry of the Arp2/3-complex. Importantly, branching 
can occur only close to the surface, as it depends on the presence of surface-bound nucleation 
promoting factors (NPF) like WASP or Act A [13j. Whether the daughter fllament becomes 
a productive member of the growing actin gel depends on the way its direction is oriented 
relative to the direction of growth and how the object is shaped. Therefore obstacle geometry 
is a crucial determinant of the resulting structural organization of the network. 



One essential function of growing actin networks is to generate force against external 
load. Force-velocity relations for growing actin networks have been measured in different 
experimental setups and for various obstacle shapes [TH - EU] . Unexpected discrepancies in 
the results indicate that obstacle shape and the resulting difference in network organization 
also has an effect on force generation. Modern electron tomography leads to visualiza- 
tion and analysis of filamentous actin networks in ever greater detail [211 E2] • While early 
electron microscopy data for the lamellipodia of fish keratocytes suggested a dendritic net- 
work of relatively short actin filaments [23], recent electron tomography data has revealed a 
more diverse structural organization, with relatively long filaments connected by few branch 
points into different and spatially extended filament subsets [22l El]. In addition, it has 
been demonstrated that the structural organization of the network strongly depends on the 
protrusion speed, with fast growing networks dominated by two symmetric diagonal filament 
orientations and slowly growing networks featuring more filaments in parallel and orthogonal 
to the leading edge [23 EH]- It is to be expected, that in the future the effect of obstacle 
geometry on network structure can be quantified by such methods as well. 

In order to achieve a complete understanding of the structure of growing actin networks, 
different modeling approaches have been developed during the last decade [27], ranging 
from microscopic ratchet models [28l |29] . rate equations for filament growth [30H32] and 
large-scale computer simulations of filament ensembles [331436] to continuum theories of 
how elastic stress propels the obstacle [3Z| and multi-scale models combining several of 
these model classes [201 I3S] . From some of these studies it has emerged that one central 
quantity characterizing the structural organization of growing actin networks is the filament 
orientation distribution, which can be directly compared with experimental results [261 ED 

MM- 

In the following, we study a model which allows us to predict the filament orientation dis- 
tribution as a function of obstacle geometry. Our reference point will be stochastic computer 
simulations of growing actin networks incorporating the molecular processes of branching, 
capping and filament polymerization [32|. For computational simplicity and deeper insight, 
these will be compared to different versions of a rate equation model [30H32] . For actin 
growth behind fiat and laterally widely spread obstacles the steady state organization has 
been predicted earlier to be either a ±35 or +70/0/— 70 degree filament orientation pattern, 
with mutually exclusive stability determined by the model parameters [311 132] • In this pa- 



per, we will extend this analysis to also predict the effect of finite size and geometry of the 
propelled object. In particular, we will derive and validate a partial differential equation, 
which is valid also for curved obstacle shapes. 

The article is organized as follows: In Sec. [IT] we introduce the model and explain how 
we analyze it. To validate our results, we compare different versions of a deterministic 
continuum model to stochastic computer simulations of filamentous networks. Subsequently, 
the impact of different piecewise straight (linear) obstacle shapes on the resulting network 



orientation patterns is analyzed in Sec. Ill, while specific examples for curved (nonlinear) 
obstacle shapes are the subject of Sec. IV Our main result is that the competition between 
two fundamentally different network architectures persists for finite-sized, piecewise linear 
and curved obstacles. 

II. MODEL DEFINITION 

Motivated by established biological observations, in our model we assume Arp2/3 to 
nucleate daughter filaments from preexisting mother filaments at a characteristic relative 
branching angle around 70° [2l [23], HI] . The exact value varies with biological system and 
analysis methods (for example, a recent value from electron tomography data for the lamel- 
lipodium of fibroblasts is 73 ± 8° [22j), but is not essential for our theory. Although the 
exact mechanism for the activation of Arp2/3-mediated branching is not yet well estab- 
lished and subject of current research, it is generally accepted that Arp2/3 is active only 
near the obstacle, where it is activated by nucleation promoting factors (NPFs) such as the 
WiskottAldrich syndrome protein (WASP) or the bacterial actin assembly-inducing protein 
(Act A) [THl 1121 113] ■ As we are interested in actin network architecture in close proximity 
to the surface of a propelled obstacle, we will focus on actin dynamics within the first few 
ten nanometers from the obstacle surface in which filament branching, capping and barbed 
end polymerization are expected to dominate actin dynamics, while filament depolymeriza- 
tion and decapping can be neglected [H]. Fig. [I|a) sketches the geometrical arrangement 
studied here. The yellow branching region within a vertical distance d^^ from the obstacle 
surface indicates the domain in which filament bound Arp2/3 is able to interact with NPFs, 
which themselves are active only close to the surface. Thus filament branching occurs in this 
branching region. Once filament barbed ends have left this domain due to their retrograde 



flow they are not able to nucleate new daughter filaments anymore and will eventually be 
outgrown by the bulk network. As those filaments subsequently do not impact the orien- 
tation distribution at the leading edge, we do not exphcitly account for their further fate 
anymore. If the lateral extension d^j. of the obstacle is large compared to the width d^^. of 
the branching region, it is possible to neglect the process of filament barbed ends growing 
out of the branching region horizontally. However, for relatively small obstacles, for instance 
viral pathogens, the lateral dimension of the propelled particle d^^ can become comparable 
to its vertical extension d^^ and the finite lateral size of the branching region has to be taken 
into account. Our main quantity of interest is the filament orientation angle 6 relative to 
the surface normal. 

In the following, we will validate our theoretical approach by comparing results from 
two complementary implementations, namely stochastic computer simulations and a rate 
equation approach, which have been introduced before for fiat obstacles of large lateral 
size [3^. Motivated by the fiat nature of the lamellipodium and also for computational 
simplicity, our modeling is restricted to two dimensions. Representative snapshots of such 
simulations for fiat and curved obstacle geometries are given in Fig.flVb) and (c), respectively. 
Briefiy, we simulate a network of infinitely stiff rods. This simplification is justified as the 
reaction kinetics that eventually determine the stable orientation distribution in steady 
state are active within a narrow branching region along the obstacle only. Because filament 
segments spanning over this region are very short compared to their persistence length 
(nanometers versus micrometers), local bending undulations are of minor importance here. 
The remainder of the filaments is considered to be entangled in the bulk actin network, 
which effectively acts as a base for the protruding filaments. However, filaments embedded 
in semidilute actin solutions have been found to exhibit substantial bending fiuctuations 
in the past [l5l |16] and therefore we effectively account for a finite uncertainty in filament 
orientation by incorporating a distribution of relative branching angles between mother and 
daughter filaments in the model as is discussed below. Each uncapped filament is growing 
deterministically at its barbed end with a fixed velocity vm- Thus we implicitly assume 
a constant density of actin monomers at the leading edge. Polymerization is quantized 
such that filaments extent by one building block of length 6f[i per unit time. Apart from 
branching, individual filaments do not interact and hence the local filament density does not 
alter polymerization. Filament barbed ends within the branching region close to the leading 



edge are possible candidates for stochastic branching and capping events. While capping is 
assumed to be a first order reaction in the number of actin barbed ends in the reaction zone, 
for branching we assume zeroth order. This is motivated by the expectation that the supply 
of activated Arp2/3 is strongly limited by the availability of NPFs at the leading edge and 
thus an effective zeroth order branching rate emerges for sufficiently high filament density 
or low capping rate |171 HH]. In the opposite limit of low filament density, the branching 
reaction is expected to become first order in the number of filaments; this effectively yields an 
autocatalytic description [30] , that is stationary only at one unique network growth velocity. 
Therefore in this limit, transitions in the filament orientation distribution are not accessible. 
However, it has been shown that both, first and zeroth order branching models, can be 
incorporated within a unified theoretical framework and that the crossover from one regime 
to the other does not have a direct impact on filament orientation in steady state [58]. As 
actin networks in most experimental setups are growing against a finite load of the obstacle 
and non-constant force- velocity curves are observed, the zeroth order branching regime is 
expected to dominate the branching kinetics eventually. 

The orientation of new filament branches is chosen from a normalized linear combina- 
tion of two Gaussian probability distributions for the angle between mother and daughter 
filament with means at ±70° and standard deviation 5°. Red filaments in the illustration 
are actively growing, while blue filaments indicate capped barbed ends that neither branch 
nor polymerize anymore and will eventually be outgrown by the bulk network and leave the 
simulation box at the bottom. In Fig. [lib) the top boundary of the yellow branching region 
defines a rigid obstacle, which excludes volume and thus prevents polymerization above its 
boundary. As a consequence the fastest filaments which are growing in close to vertical 
direction are stalled by the obstacle. The vertical velocity of the obstacle is a parameter of 
the model and set to a constant value t>nw directly. This growth velocity should be regarded 
as the effective outcome in steady state, determined by the details of filament-obstacle inter- 
actions (including regulation by biochemical factors). Steady state growth is possible within 
the velocity range, < tinw < ^fii- Under these conditions a well defined filament number 
and orientation distribution evolves. 

As a powerful analytical alternative to the stochastic framework introduced above, we 
also develop a deterministic rate equation for dendritic network growth based on earlier 
approaches of this kind PDH52] . The evolution of the orientation distribution N{6,t) of 



uncapped filament barbed ends integrated over the whole branching region evolves in time 
by branching and capping events and by filaments growing out of the branching region. This 
translates into the following ordinary differential equation: 

^^^^ = h [ W(e, 9')N{9', t) d9' - hNi9, t) - k^,{9, v^^)N{9, t) . (1) 

Ot J^ ^ V ' V ' 

^^ V ' capping outgrowth 

branching 

Here again, capping is assumed to be a first order reaction proportional to the number of 
existing filaments with the proportionality constant, a reaction rate kc. The probability 
of branching at a given angle is determined by the weighting factor distribution, yV{9,9'), 
which is modeled as a linear combination of Gaussians with maxima at branching angles 
±70° between filaments and standard deviation cr = 5°. ^ is the filament orientation angle 
measured relative to the vertical direction. The branching reaction is independent of the 
absolute number of existing filament barbed ends in this model (i.e. a zeroth order reaction). 
Hence, the branching rate k^ is normalized by the total number of new filament ends, 

^b = ^ , >Vtot = / / W{9, 9')N{9', t) A9' d9 . (2) 

I^Vtot J J 

The reaction rate fcb indicates the number of new branches per unit time. Most importantly 
for the specific scope of this work, obstacle velocity and shape enter Eq. ([I]) via the outgrowth 
rate fcgr. As filaments are polymerizing with a constant velocity fgi in their individual 
direction, some filaments will not be able to keep up with the moving obstacle and thus 
leave the branching region. The precise expression for the rate of outgrowth fcgr strongly 
depends on obstacle geometry and alters the resulting steady state orientation distributions 
as will be analyzed in detail in the following. 

III. LINEAR OBSTACLE SHAPE 

We begin our analysis with linear obstacle shapes as illustrated in Fig. [l|a) and (b). The 
obstacle (and therefore the leading edge of the network) moves with a constant velocity 
fnw towards the top. Our region of interest, the branching region, extends vertically to a 
distance d^^ from the leading edge of the network. The lateral obstacle width is indicated 
by d^^. 



A. Flat large obstacle 

For obstacles that are laterally widely extended, for instance the leading edge at the front 
of the lamellipodium in a migrating cell, we have d^^ ^ d^^. and thus horizontal filament 
outgrowth can be neglected locally. Mathematically, this means that we can use periodic 
boundary conditions in the horizontal direction. This case has been analyzed before in a 
similar manner as done below for finite-sized and curved obstacles [32], and therefore we 
recapitulate the most important results as a reference case. In this simple case, outgrowth 
of filaments from the branching region can only occur in negative normal direction. While 
the network grows with velocity v^^ in positive normal direction, single filaments grow with 
velocity f m in their individual direction. The projected polymerization velocity depends on 
the filament orientation, as 

^m(^) = ^^fiicos6'. (3) 

Filaments having a larger absolute orientation than the critical angle 

9c = arccos I -^ ] , (4) 

are thus not able to keep up with the speed of the obstacle and are subject to outgrowth 
with a rate 

{0 if 1^1 < 9c 

(5) 
"-(g%7' if|^|>^c 

The factor 2 results from the assumption that new filaments branch-off from existing fila- 
ments on average in the center of the branching region of vertical width d^^.. Outgrowth 
vanishes for filaments with an orientation smaller than a critical angle 9c and from there 
increases to its maximum at \9\ = 180°. 

We solve Eq. ([I| numerically by introducing 360 angle bins and iterating the equations 
until a steady state is achieved, which then can be compared to the results of the stochastic 
computer simulations. In order to achieve a deeper understanding, we also introduced a 
coarse-grained version of Eq. ([I]) that can be treated analytically. In Eq. (fTl), the number of 
filament ends in the branching region with angles between 9 and 9 + d9 is given by N{9, t)d9. 
By extending this integration to sufficiently large (A^ = 35°) angle bins 

e+A6»/2 

Ne= J N{9\t)d9', (6) 

e-A6»/2 



and further assuming that branching is restricted to pairs of angle bins with a relative angle 
difference of 70°, a system of five coupled ordinary differential equations results: 

aiV_70o 1 : 



hNoo - {K + k^A70°)) iV„7oo (7) 

hN+35^ - {h + k^A35°)) iV_35o (8) 

k{N.70o+N+jo°)-k,Noo (9) 

fcbiV-35° - {kc + k^^{35°)) iV+35o (10) 

hNoo - (fee + fc^,(70°)) iV+7oo , (11) 



>9iV-35° _ 1 

dt 2 

dNoo _ 1 

dt ~ 2 

dN+s,o ^ 1 

dt 2 

(9A^+7o° 1 



with 



kh 



dt 

fcb k 



>Vtot Ar_7oo + Ar_35o + Noo + Ar+350 + Ar+7o° 
Here we have also assumed that branching of filaments with orientations \6\ > 87.5° can 
be neglected. The five equations are symmetric around 0°, and thus only three of them 
are independent. Nonlinearity and the coupling of all five equations is introduced by the 
branching term due to the zeroth order branching reaction. 



By algebraically solving Eq. rt7|)-(Eq. (11) for the stationary state, we obtain two phys- 



ically meaningful solutions. The first solution, 

NH^Io = 

j\Tss35 _ u 1 

JV_35o - «'b4(fc^_|.fc±(35o)) 

ATq-ss = (12) 

/yssSS u 1 

^^+35° '^b4(fc^+fc±(35o)) 

iv;f5o = 0, 

represents a dominant ±35 degrees orientation distribution in the steady state while the 
second solution, 

70 _ , fcc+fcg^.(70°)-^2fc,(fc,+fc^^,(70°)) 
^ -™° ~ ""^ 2(fe,^,^70°)-fe?) 

NH™ = 
JVfssYO _ /. / V 2fcc (13) 



,.ss70 _ , fcc+fc^,(70°)-y2fcc(fcc+fcgL(70°)) 
iV+7oo - fcb 2(fcJ,2(70°)-fc2) 
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corresponds to a competing +70/0/— 70 pattern. 

The stability of these two fixed points can be analyzed with respect to the parameters 
fcb, kc, ^^(35°), fc^(70°) and d^^ using linear stability analysis. The result of this analysis is 
independent of the branching rate fcb and therefore this parameter has no infiuence on the 



stability of the system. According to Eq. (12) and Eq. (13), however, the total number of 
filaments in steady state is proportional to A;b and therefore this parameter can be used to 
adjust the model to experimentally measured densities. The two parameters A;j"j,(35°) and 
/c^(70°) are not independent, but rather both of them are determined by the obstacle velocity 
fnw as given in Eq. (p|. In the following, we will omit the ill-defined cases k^ = fc^(35°) = 
and kc = fci,(70°) = 0, for which the filament number diverges. We find that the stability 
of both fixed points changes when 



fccfcj-r(70°) = A;2 + 4fcefcj-r(35°) + 2fcj-/(35° 



(14) 



and that either one is asymptotically stable, while the other is a saddle. Thus the simple 
model suggests that exactly two possible network architectures exist with mutually exclusive 
stability. 



Eq. ( 14 ) can now be used to obtain the respective network velocity Vnw of the transition. 



For a critical angle 9c > 70° (i.e. k^j.{70°) = ^^(35°) = 0), Eq. (14) is never satisfied (note 
that kc > 0). For 35° < ^c < 70° (i.e. ^^^(70°) > A A;^,(35°) = 0), 



Vn^ = ^Kdi + Vf,i cos(70°) , for 35° < ^e < 70° 



(15) 



satisfies Eq. (14). For network velocities with a critical angle 6c < 35° (i.e. fci.(70°) > 
A;j^(35°) > 0), two solutions emerge. 



^^nwi,2 = I (-3/i;c<i^r + 8^fiicos(35°)) 



±|^Mbr (^cO?br + 16vfii cos(35°) - 16wfii cos(70°)) 
for Or < 35° . 



(16) 



It turns out that solution Eq. (15) is valid for network bulk velocities Wfiicos(35°) > fnw > 



ffiicos(70°), while for solution Eq. (16) to be valid fnw > fmcos(35°) has to be satisfied. 



This is never the case for the negative square-root in Eq. (16) and so we can neglect this 
solution in the following. In order to further simplify our equations, we define the reference 



velocity Uc = kcd^^ resulting from the capping rate. Eq. (15) and Eq. (16) then become 

■y^w 1 -u, 



^^fii 
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+ cos(70°) 



(17) 



and, 

i(-3g^ + 8cos(35°: 



''HI 



+ lJ^. (^ + 16cos(35°) - 16cos(70°; 



(18) 



respectively. Thus the capping rate emerges as an essential parameter through the effective 
velocity Uc- 

In Fig. p[a) the stability diagram from the analytical analysis is illustrated, showing 
the regions in which either the ±35 degree distribution is asymptotically stable and the 
+70/0/— 70 pattern is a saddle or vice versa. In order to compare these results to the full 
numerical solution of Eq. ([I]), we define the relative difference of filaments in the angle bin 
around 0° to the one around 35° as an appropriate order parameter in steady state. 

For a perfect +70/0/— 70 distribution this parameter will approach +1, while for the com- 
peting ±35 pattern it will approach —1. The transition between the two patterns is defined 
when the order parameter changes its sign. According to this definition of the transition 
point, it is now possible to numerically solve the continuum model Eq. ([I| and to classify 
each observed stationary state as one of the two phases. These results are presented as a 
phase diagram in Fig. W[d) and were obtained for different branching rates /cb, which indeed 
has no significant influence as expected from the analytical considerations. In the solution 
of the full continuum model the capping rate only has a very limited influence on the stripe- 
like pattern of the phase diagram. Contrarily, in the phase diagram resulting from linear 
stability analysis (Fig. [2|^a)), the ±35 pattern vanishes at large capping rate. This artifact 
of the reduced rate equation model is introduced by the assumption that filaments with an 
orientation larger than 87.5° do not branch. As the full rate equation does not share this 
assumption, it does not predict the elimination of the ±35 pattern for large capping rate. 

B. Flat finite-sized obstacle 

For flat, but finite-sized obstacle shape, filaments might also leave the branching region 
horizontally. Therefore a second outgrowth rate needs to be incorporated into the rate 
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equation (Eq. ([T])): 

{ f nw tan 6> -r \n\ ^ n 
(C./2) ' ' - \ (20) 

(d=/2) 11 l^^l > (^c 

We find that for tliis scenario, it is still possible to analyze the stability of the fixed points 
of Eq. ([7])-Eq. (11). The only difference that accounts for the change in obstacle geometry 



at this point is that a combination of outgrowth rates k^riO) = k^ri^) + |^^(^)| replaces 
the term for exclusively orthogonal outgrowth /ci(6') considered in the previous section. For 
filaments growing in the A^o° angle bin, both outgrowth rates vanish, because these filaments 
are growing at a high enough velocity parallel to the lateral boundaries of the obstacle. 
The same arguments as used before also apply for the stability analysis of the steady 



states here. Hence, we can begin with Eq. ( 14 ) and evaluate the stability of the two steady 



states according to the adjusted outgrowth rate. To determine the transitions, we will follow 
a similar strategy as before. Starting from small network velocities v^wj we will treat the 
different possibilities for the outgrowth rates given in Eq. ([s]) and Eq. (20) in a case-by-case 
analysis. 

For 6c > 70°, the orthogonal outgrowth rates vanish, as all relevant filament orientations 
grow faster towards the top than the obstacle and are slowed down to fnw in this direction. 
The outgrowth rates for the different orientation bins therefore read. 



f „w tan 61 



k^M = and k={e) = "" /^x for 6 = [35°, 70°] . (21) 



Inserting kg,- = fcr, + k in Eq. ( 14) and omitting again the ill-defined cases of /cb = 0, fee = 



and d^j. = here (and additionally d^^. = in the following), it can be shown that Eq. (14) 
is never satisfied. 

Increasing the network velocity to the point where 35° < 6c < 70° leads to outgrowth 

rates, 



t^fii cos (70 

in orthogonal direction according to Eq. ^ and, 

t;, ^tan(35°) _ t^fiisin(70 ° 



ki^i35°) = and ^^^,(70°) = °" .,' ^ ^ , (22) 



;t;(35°) = ^-;rr; ^ and fc;(70°) = '''']Z)n. ' (23) 



in the lateral direction according to Eq. ( 20 ) . These rates yield a quadratic equation for the 
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network velocities f nw that satisfy Eq. ( 14 ) 



;tan2(35°)^ 



8tan(35°)^ 2^ 



+ 



1 - 2t;fiisin(70°)^ + 2t;fiicos(70°)^ 



(24) 



for 35° < ^c < 70° . 



Here, the reference velocity, Uc = kcd^,., together with a new parameter for the length scale 
ratio in lateral and orthogonal direction, rbr = '^br/'^br' ^^^ identified to determine stability. 
For a critical angle 6c < 35°, the adjusted outgrowth rates read. 



kiA0) = '"^";;T'^ and kZiO) = ^S^ for 6 = [35°, 70°] 



(dim 



(rfb=/2) 



(25) 



These rates together with Eq. ( 14 ) and the effective parameters Uc and rbr can be simplified 



to a second quadratic equation defining the network velocities at the transitions. 



+ Vn 



6^ - 16t;mcos(35°)^ + 16t;fiisin(35°)^ 



+ 



1 + 2t;fiicos(70°)^ - 2i;fiisin(70°)^ - 8t;mcos(35°)^ + 8t;fiisin(35°)^ 



+ 8<cos2(35°)^ - 16<cos(35°) sin(35°)^ + 8<sin2(35°)^ 
for Or < 35° . 



(26) 



According to Eq. (24) and Eq. (26), the regions of stability for the two different stationary 

As the finite obstacle 



orientation patterns in parameter space are illustrated in Fig. p, 

width introduces an additional independent parameter rbr, the full parameter space is now 



three dimensional. Network growth velocities v^w that fulfill Eq. (24) together with the 



condition t>fiicos(35°) > v^^ > t'fiicos(70°) or Eq. (26) at v^^^ > ffiicos(35°) are identified 
as transition points between the different orientation distributions in the parameter space 
spanned by kc, rbr and fnw In the limit of large length scale ratio, rbr — ^ oo (i.e. d^^ 3> d^^), 
lateral outgrowth can be neglected and the phase diagram approaches the results for periodic 
boundary conditions as given in Fig. [2Fa). For relatively small rbr ^ 7, the lateral outgrowth 
of intermediate filament orientations at around ±35° is sufficiently large to prevent this 
orientation pattern from being stable in this model, independent of the growth velocity (cf. 
Fig.gc)). 



Using the order parameter Eq. (19), we numerically sampled the parameter space accord- 



ing to Eq. ([I]) with adjusted outgrowth. The isosurface C(fnw, kcr^j-) = is extracted from 
the three dimensional data and shown in Fig. W[e). In the limit of large length scale ratio, the 
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results coincide well with the case of periodic boundary conditions (cf. Fig. [2|^d)). It is also 
confirmed that for small values of rtr the ±35 pattern is not stable anymore at intermediate 
velocities (cf. Fig. |2Vf)). Thus the overall agreement between the simple analytical model 
and the full numerical solution of the rate equation approach is surprisingly good. 

Fig. [3] compares filament orientation distributions obtained in steady state from stochastic 
network simulations, the full continuum model and the simplified continuum model. The 
parameters of the two cases shown in (a) and (b) do not differ in their network velocity v^^, 
but rather only in the obstacle geometry, in this case given by length scale ratio rbr = d^^/d^^. 
For small rbr = 3 at intermediate velocities, the ±35 pattern of the network is not stable 
anymore and rather the network organizes in the alternative ±70/0/— 70 distribution (cf. 
Fig. [3ta)). For larger ratio rbr = 20, the network organizes in the ±35 degree distribution (cf. 
Fig. |3](b)). In the stochastic simulations, outgrowth rates are not explicitly incorporated, 
but rather emerge as a direct consequence of the obstacle geometry. Thus the computer 
simulations nicely validate both, the full and reduced continuum approaches, and all three 
seem to capture the essential physical mechanisms determining network structure. 

C. Obstacle with tilted straight contour 

It is not trivial to find an explicit expression for the outgrowth rate out of the branching 
region behind a skewed linear obstacle, which is rotated according to a constant skew angle 
(y9, with —90° < if < ±90°. However, for reasonably small skew angle ip and width of 
the orthogonal branching region d^^ an approximation can be obtained within a rotated 
coordinate frame, that is adjusted to the constant obstacle skew. Fig. |4Fa) shows a sketch 
of the obstacle geometry we are interested in at this point. The branching region is given 
by a rhomboid, with rotated (or skewed) upper and lower side. The orthogonal width of 
the branching region d^^ is again defined parallel to the lateral sides of the obstacle, while 
its horizontal width is given by d"^^. Network velocity v^w is defined as before in the vertical 
direction. To find an explicit expression for the rates of filament outgrowth, a rotation of 
the coordinate frame to the point where the skew angle yj of the branching region vanishes 
simplifies the situation (Fig. |4](b)). In this frame the obstacle appears fiat horizontally, very 



similar to the setup analyzed before in Sec. IIIB The finite skew angle ip, however, manifests 



in a non- vanishing lateral network growth velocity f„^, while the vertical propulsion speed 



15 



of the obstacle is the orthogonal part of the network velocity in this frame, v:^^. (Where 
applicable, we will denote variables defined in the rotated coordinate frame by ~.) In the 
following we will refer to the original coordinate system Fig. |4](a) as the lab frame and to 
the rotated picture Fig. mh) as the obstacle frame. 

Using the transformation to the obstacle frame for a given skew angle ip, we can now 
deduce the relevant parameters. The filament orientation angles 6 are given relative to the 
obstacle, i.e. for 6 = 0° filaments are growing vertically in the rotated frame. The two 
components of obstacle growth velocity are given by, 

^nw(V') = ^'nw COS if and v=^{'f) = Wnw slu (f . (27) 

The two obstacle diameters, which scale the vertical and horizontal outgrowth rates are, 

du ((p) = dLcosif and d^ (ip) = — ^ . (28) 

COS(p 

According to these definitions, we can now reformulate the continuum description in the 
obstacle frame. The rate equation model Eq. ([I]) for the temporal evolution of orientation 
dependent filament density remains unchanged with respect to the transformed orientation 
angles 9, however with adjusted outgrowth rate fcgr(6',fnw) given by the sum of orthogonal 
outgrowth, 

~, ^ [o if 1^1 <^c 

k^,{e,^,Vr,^)={ , , (29) 

[ (J4(^)/2) II |f I > fc 

and the absolute value of lateral outgrowth. 



"nw(y)+^ifw(y)tane 



if 1^1 < 0, 



KA0,^,Vr..) = l (^-(^)/^l ''- , (30) 



"nw(y)+^flisin6> 



if 1^1 > e. 



with. 



arccos ( "^^^ ] . (31) 



^^fii 



The lateral outgrowth rate compared to Eq. (20) is biased by the finite horizontal velocity of 



the network in the obstacle frame wj^. This additional feature also breaks the symmetry in 
the resulting steady state filament orientation distributions. The sign of the outgrowth rate 
determines, whether filaments are growing out to the left (fc^ < 0) or the right (fc^ > 0) 
of the branching region, which is only relevant when keeping track of the filaments' fate 
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after outgrowth. This will become important later in Sec. IV A but at this point the sign 
of lateral outgrowth is suppressed. 

To test our approximation, we again compared steady state filament orientation distribu- 
tions, obtained by numerical iteration of Eq. dl]) within the obstacle frame and subsequent 
transformation to the lab frame, against stochastic simulations. Typical cases are shown in 
Fig. [5Fa)-(d) for a moderate obstacle skew angle oi ip = 20°, fast and slow network velocity 
in the lab frame fnw? and small and large horizontal obstacle diameter as indicated by the 
length scale ratio rbr (Fig. IS] also shows results obtained with a PDE-model which we will 
introduce below). If the results are interpreted directly in the adjusted obstacle frame, very 
similar orientation patterns as the familiar ±35 and +70/0/— 70 peaked distributions are 
realized in steady state. In the limit of very large length scale ratio rbr, outgrowth in the 
lateral direction can be neglected again and the resulting patterns resemble orientation pat- 
terns for a network growing behind a large obstacle, only now within the rotated obstacle 
frame (cf. Fig. [stb) & (d)). For small rbr in combination with a relatively high obstacle 
velocity, horizontal outgrowth cannot be neglected and the orientation distribution has to 
be interpreted in the lab frame (Fig. [sla)). In between a transition occurs where the solu- 
tions from stochastic simulations already have to be interpreted in the obstacle frame, while 
results from the rate equation model are still to be interpreted directly in the lab frame 
(Fig. |5](c)). The reason for this is the spatial resolution in the lateral filament position, 
which is introduced in filament based stochastic simulations and neglected in the continuum 
model. Laterally nonuniform spatial filament distributions lead to differences in the effective 
outgrowth of filaments to the sides of the branching region. In order to include a similar 
spatial resolution in the continuum model, in the following we will incorporate an additional 
advection term in the equation. In this way it will also become possible to treat nonlinear 
obstacle shapes within the continuum description. 

IV. NONLINEAR OBSTACLE SHAPE 

In this section, we focus on actin network dynamics in the tail of curved (i.e. nonlinear) 
obstacle geometries. The definitions of the relevant parameters stay the same in this geome- 
try as introduced in Fig. U[a). Fig. m c) features a representative steady state network from 
stochastic simulation. In the following, we will stepwise extend the continuum model Eq. ([I]) 
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with the goal to also incorporate nonlinear obstacle shapes in the description. We will do 
this first within an adjusted ordinary differential equation (ODE) model in a piecewise-linear 
approximation of the obstacle surface and then within a partial differential equation (PDE) 
model, where an advection term will explicitly govern lateral filament growth in the rate 
equation, additional to the reaction terms. 

A. Piecewise-linear obstacle approximation 

Nonlinear obstacle shapes can be treated within a piecewise-linear approximation for the 
obstacle surface in which the ODE-model from above is still applicable. The construction of 
this approximation is sketched in Fig. |4Fd)-(f) for a parabola-like curved obstacle surface. 
The branching region is divided laterally in n sections of equal size d^^^ = d^^/n {i = 
1 , . . . , n) . The orthogonal width of the branching region in each section remains unchanged 
at d^^- = d^^ . The resulting obstacle approximation is a combination of lateral sections 
with skewed linear shape according to skew angles ipi, very similar to the obstacles discussed 
before. The individual sections (i.e. the ODE-model equation of each section) are coupled 
to their direct neighbors via the lateral outgrowth rates fc^ ^ to the left and right boundaries 
as indicated in the sketch. In case of periodic lateral boundary conditions of the obstacle, 
sections i = 1 and i = n are also coupled via outgrowth. Individual outgrowth rates within 



the local obstacle frame of each section can then be formulated according to Eq. (29) and 



Eq. (30) as before. 

Assuming a single linear obstacle is the simplest approximation of a nonlinear geometry 
(Fig. pFd)). To increase accuracy, the obstacle can be subdivided laterally in ever smaller 
subsections. Due to the left-right symmetry of many biologically relevant obstacle shapes 
(e.g. the parabola-like shape in the sketch) it is sufficient to obtain results in one half- 
space of the obstacle geometry. This translates into a problem of n/2 coupled ODEs for an 
even number of subsections in the piecewise-linear approximation. In order to exploit this 
symmetry, lateral outgrowth from the section adjacent to the middle of the box has to be 
adjusted: Filaments that would leave the section in the middle towards the other side of 
the obstacle are reinjected again at this position with mirrored orientation, 6 -^ —9 . For 
instance in case of a triangular approximation of the nonlinear obstacle (Fig. pt^)), only 
the outgrowth behavior at one boundary needs adjustment compared to the problem of a 
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skewed linear obstacle shape that was already discussed before. Fig. |6]^a) and Fig. [TJ^a), 
show steady state filament orientation distributions for a laterally nonperiodic parabola and 
periodic-cosine shape obtained by the ODE-model in triangular approximation respectively. 
The nonlinear obstacle shapes together with their respective triangular approximations are 
illustrated in Fig. W[c). To parametrize nonlinear parabola and cosine geometries in this 
context, we are using the left hand side skew angle of the triangular obstacle approximation 
if in combination with the horizontal width of the obstacle indicated by rbr- 

B. Reaction-advection equation 

The continuum limit of an infinitely large number of piecewise-linear subsections, each 
of infinitesimal width, yields a PDE-model for nonlinear obstacle shapes. The model equa- 
tions in this limit are complemented by an additional filament advection term, similar to 
hydrodynamic models, that incorporates lateral filament growth in Eq. ([I]), 



dt dx J (32) 

-hN{9, X, t) - k^,{6, if, v^^)N{9, X, t) , 

with the horizontal filament growth speed, v'^^iO.Lp.v^^), and the vertical outgrowth rate in 
the lab frame, ^^(6*, v?, fnw)- Both, filament growth velocity and outgrowth rate are functions 
of filament orientation 9 as well as the local obstacle skew angle ^{x), which is active at 
lateral position x. Here again, we average outgrowth of filaments in the vertical direction 
by an effective rate, while horizontal growth is now spatially resolved due to the additional 
advection term. 

For arbitrary obstacle shape, that can be expressed in terms of an analytic function o(x), 
the local skew angle can be written as, 

, , f doix)\ ,„„, 

(p[x) = — arctan I — — I . (33 j 

The two components of the filament growth velocities within the local obstacle frame can 
thus be written as, 

if 1^1 < ^c 

^m(^,<^,t'nw) = { „ . . ' (34) 

^nw(v^) - ^fii COS 9 if 16*1 > 9c 
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and, 



^filV^5 



V5,^'i, 



v=-„i^)+vmsm9 if |^| > ^c 



(35) 



with the components of the obstacle velocity, v^^{ip) and v^^{(p), defined as in Eq. (27) and 



6c as in Eq. (31). A transformation (i.e. rotation) from the local obstacle frame into the 



local lab frame at each horizontal position along the obstacle surface subsequently leads to 



filament growth velocities, which enter the PDE-model Eq. (32) due to vertical outgrowth 
and advection, 



(36) 



V^^{6,(p,Vnw) = V^^{6, (f, Vnw) COS !f + Vfi^{e,ip,Vnw) simp 

Vfil(6',V2, fnw) = - V^i{9,ip,Vnw) simp + V^^{9,(p,Vn^) COS ip. 

Note, that filaments with a positive v^i{9, p) are growing vertically in opposite direction to 



the obstacle. Using these expressions, we can now rewrite the advection term in Eq. (32) 
as. 



d{N -n 



m) 



.ON 



dv 



dx 



^^ dx dx 



fii 



^=dJl^^dv^dp 
^^ dx dp> dx 



with, 



dv, 



fii 



dp 



d{viisinp) d{v^^cosp) 



dp 



dp 



(37) 
(38) 



and substituting Eq. (34), 
d (-Cgi sin p) I 



if 101 < Or 



dp 



— fnw sin^ p + Vnvj cos^ p — Vfn COS 9 COS p — vm sin psiia9 if |6'| > 6'c 



(39) 



and Eq. (35) 



d{v, 



'fll cos Ifij 



dip 

f nw cos^ p — Wnw siu^ p — 2v-aw COS p siu p tau 9 — f nw cos^ p cos~^ 9 if 1 6*1 < 6'c (40) 

Wnw cos^ p — f nw siu^ p — v^\ siu 9s\n.p — v^i COS p COS 9 if 1^1 > 6'c 



and using Eq. ( 33 ) , 



dp 
dx 



d^ojx) 



l+(^ 



(41) 



As we average outgrowth in the vertical direction by the rate k-^-^{9, p, v^w) in the lab frame, 
while horizontal growth is now spatially resolved, we have to take into account the finite lat- 
eral movement of filaments along the obstacle geometry, when considering filaments growing 
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out of the branching region vertically. To approximate the filament outgrowth rate accord- 
ingly, we again assume a locally linear obstacle shape at each horizontal position to find an 



analytical expression for the outgrowth rate in Eq. (|32|) 



\±^ , . (42) 



To calculate steady state filament orientation distributions from this PDE-model, 



Eq. (32)-Eq. (42) can be numerically iterated in discretized space and time using for in- 
stance a second order upwind scheme for the spatial differential in combination with an 
Euler method for temporal iteration [l9]. For this procedure, we are using a uniform fila- 
ment distribution in space and orientation as initial condition and, dependent on the specific 
obstacle shape considered, either periodic boundary conditions or zero filament density at 
the respective infiow boundary. 

C. Linear obstacle shape revisited 



The important advantage of using the PDE-model Eq. (32) compared to the initial ODE 
Eq. ([I]) to solve for steady state filament patterns clearly lies in its additional spatial reso- 
lution. This benefit not only increases the applicability of the model to a broader range of 
obstacle shapes, but also manifests itself when solving for specific parameter combinations 
for simple nonperiodic linear fiat obstacles that have been already accessible using the ini- 
tial ODE equations. Fig. [ste) illustrates the steady state filament orientation distributions 
averaged over the whole branching region behind an obstacle of small width (i.e. rbr = 3). 
For such laterally small and fiat obstacles, the ODE-model predicted the absence of the ±35 
orientation pattern for arbitrary obstacle velocities v^-w (cf. Fig. [2]). Due to their additional 
spatial resolution horizontally, the PDE-model as well as results from stochastic network 
simulations yield this orientation distribution nevertheless for a small range of parameters 
at relatively slow obstacle velocity, fnw/'^fii = 0.2. As shown in Fig. ^i), the spatial steady 
state filament distributions from the PDE-model are nonuniform horizontally and in good 
agreement to results from stochastic simulations. This has an impact on the effective lateral 
outgrowth of filaments in the model and thus is able to change the final (spatially averaged) 
filament distributions characteristically as shown in this example. 
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D. Parabolic and periodic-cosine obstacle shapes 

The PDE-model allows us to systematically analyze actin networks in the tail of nonlinear 
obstacle shapes. As a proof of principle, in the following we have chosen two different obstacle 
shapes which are motivated by highly relevant biological examples. On the one hand, we 
model a laterally nonperiodic parabolic obstacle shape (blue solid line in Fig. |4](c)) which 
is a first approximation for the spherical or ellipsoidal shape of actin-propelled intracellular 
pathogens. On the other hand, actin growth behind a periodic-cosine geometry (red solid 
line in Fig. |4](c)) is analyzed, motivated by the laterally widely spread but corrugated leading 
edge of a crawling or spreading cell. 

In Fig. [6] and Fig. [7| the resulting steady state patterns from the three different model- 
ing approaches (PDE-model, ODE-model in triangular piecewise-linear approximation and 
stochastic network simulations) are shown for specific parameter combinations. In (a) the 
filament orientation distributions spatially integrated over the left side of the symmetric 
obstacle are shown. Here, the PDE-model corresponds very well to stochastic network simu- 
lations. Despite the rough piecewise-linear approximation of the nonlinear obstacle shape in 
the ODE-model, the results capture the overall characteristic of the resulting network pattern 
very well. Subfigures (b) illustrate the resulting spatial distributions of certain characteristic 
filament angles. For specific parameter cases separate spatial domains emerge over the hori- 
zontal obstacle dimension, in which the signature of the two alternative steady state network 
patterns, i.e. +70/0/— 70 and ±35, are observed. For a better overview over the complete 
spatially resolved orientation patterns from the PDE-model, Fig. ^c) and Fig. m c) illustrate 
the two dimensional filament distribution in heat plots, where regions of cooler color indi- 
cate higher filament density. The apparent coexistence of alternative patterns resolves when 
additional dotted lines are included in the heat map to indicate local filament orientations 
that correspond to the characteristic orientations ^ = 0°, ±35°, ±70° within the obstacle 
frame, locally orthogonal to the obstacle surface along its contour. Although coexistence of 
alternative filament orientations is present in the lab frame, an interpretation in the local 
obstacle frame yields one unique pattern only that is present along the lateral extension of 
the obstacle. 
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V. DISCUSSION 

In this article, we have introduced several theoretical approaches at different levels of 
complexity with the common goal to predict actin filament orientation patterns at the sur- 
face of stiff two dimensional obstacles, whose surfaces promote growth of actin networks. 
This problem is not only central for many important health and disease related biological 
phenomena, such as migration of animal cells or propagation of pathogens in their host, it 
is also a prominent example for a physical process whose properties emerge on a mesoscopic 
length scale in relative independence of the details of the underlying molecular processes. 
For example, the competition of the two fundamentally different network architectures dis- 
cussed here does not rely on the exact value of the branching angle, but rather on the fact 
that any branching angle below 90 degrees can lead to the possibility that fundamentally 
different network architectures satisfy the simultaneous requirements of forward growth and 
side branching. Our model focuses on the geometrical aspects of this situation (in particu- 
lar filament orientation and obstacle shape). In the future, it might be extended by other 
important aspects of this complex biological system, including the details of the filament- 
membrane interaction and the role of filament bending [201 138] . 

In general, all our results were benchmarked against stochastic computer simulations 
which in principle can be used to include many details of the underlying molecular processes. 
Here we have adopted the established view that polymerization, branching and capping are 
the dominating processes in the context of growing actin gels. For relatively simple linear 
obstacle shapes, a rate equation (ODE) model yielded accurate results and, due to its reduced 
complexity, also allowed for analytical progress. Using this approach, we found that either of 
two competing orientation patterns emerges in steady state and transitions between the two 
are triggered mostly by changes in the velocity of the obstacle. A similar result has already 
been obtained in earlier studies for a fiat obstacle with large lateral extension, where filament 
orientation distributions with characteristic peaks at either +70/0/— 70 or ±35 have been 
found to dominate the steady state of growing networks [3ll [32] . These results were now 
extended to finite sized obstacles. We find that for very small objects, the +70/0/— 70 pattern 
is dominant, while for larger objects, mutual stability of the two patterns recovers. In case 
of obstacles with a tilted straight contour (linear obstacles), again very similar patterns are 
predicted. However, for most parameter configurations they have to be interpreted within 
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a rotated obstacle frame. 

Curved (nonlinear) obstacle shapes have been analyzed in the ODE-model using a 
piecewise-linear approximation of the given geometry. As an alternative, a continuous spatial 
coordinate in the horizontal direction was introduced to the model equation by incorporat- 
ing a filament advection term, thus yielding a PDE-model. The additional benefit of this 
spatial resolution clearly lies in the gain of accuracy compared to results obtained earlier. 
For specific parameter combinations, it was shown that the lack of spatial resolution can 
lead to incorrect predictions regarding the resulting filament patterns in steady state. Anal- 
ysis of nonlinear obstacle shapes in the PDE-model yielded an apparent coexistence of the 
two competing network patterns at different horizontal positions along the obstacle surface 
in the lab frame. A transformation to the obstacle frame, locally at each lateral position, 
however revealed that one unique orientation pattern is dominating the network structure 
as before. Again, these results are in good agreement with computer simulations. 

In order to experimentally test these theoretical predictions, the method of choice would 
be electron tomography [221 1211 [2S], although in the future super-resolution microscopy 
like dual-objective STORM might become an interesting alternative [50]. In the context 
of a rapidly increasing image quality of electron microscopy (EM) data for actin networks, 
the correct analysis and interpretation of the measured observables becomes increasingly 
important. Therefore, a detailed understanding of the structural network characteristics 
emerging under different situations is indispensable. In this work, we have also shown 
that the extracted orientation distribution of actin filament networks from experimental 
EM images growing at the surface of nonuniform obstacles can be easily misinterpreted. 
We have shown that a standard measurement in the lab frame would yield apparently novel 
patterns that could not be explained by existing models. For a correct interpretation of such 
findings, we introduced a rotated obstacle frame which lies locally orthogonal to the obstacle 
surface. Within this reference frame, the seemingly novel filament orientation patterns are 
rationalized and can be understood in terms of the two well known filament distributions, 
peaked at +70/0/— 70 or ±35. Thus our results also contribute to improving the way 
experimental data can be compared to theoretical predictions. 

A convenient setup to test our predictions would be the combination of electron to- 
mography with biomimetic assays with particles of various shapes, which earlier have been 
analyzed mainly in regard to macroscopic variables such as propulsion velocity and shape of 
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the comet tail [TT]. Apart from obstacle shape, the obstacle velocity relative to the filament 
polymerization speed is a central parameter in the model which triggers transitions between 
filament orientation patterns. This velocity could be adjusted for instance by applying a 
force against network growth diminishing steady state growth. Changing the biochemical 
reaction rates for capping and branching is expected to be of minor importance and could 
be tested by varying the concentration of purified protein solutions in the assay. As an addi- 
tional benefit of introducing different obstacle geometries in experiment, this might also lead 
to a deeper understanding of Arp2/3 activation close to the surface of the obstacle, which 
is one of the most important questions still to be clarified in the context of actin-driven 
motility. 
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FIG. 1: Model of a growing actin network behind a rigid obstacle, (a) Sketch of the setup. 
Branching can occur only in the yellow region. The main quantity of interest is the filament 
orientation angle 9 relative to the surface normal, (b) Snapshot of a stochastic simulation based on 
the reactions of individual filaments for straight (linear) obstacle shape, (c) Snapshot of network 
growth behind a curved (nonlinear) obstacle shape. The network is growing in two dimensions. 
Red filaments are actively growing at their barbed end in direction 9. Blue filaments have been 
capped and will eventually be outgrown by the network and leave the box at the bottom. Growth 
of the fastest filaments is stalled by the obstacle. 
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FIG. 2: Phase diagram for flat finite-sized obstacles, (a-c) Analytical results from linear stability 
analysis of the simplified continuum model with 5 angle bins, (d-f) Full numerical solution of the 
rate equation model with 360 angle bins. In (a) and (d), 2D slices through the 3D hypersurface are 
plotted within the kc-Vnw plane in the limit of rbr -^ oo, which corresponds to periodic boundary 
conditions or very large objects, (b) and (e) show the hyperplane of transition points between the 
two filament patterns in the full 3D parameter space. In (c) and (f), slice plots in the rbr-i'nw plane 
are shown along the transparent gray plane in (b) and (e) at kc = 0.05. 
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FIG. 3: Comparison of the steady state orientation distributions for actin networks growing in the 
tail of a flat obstacle of finite size at fnw/f^m = 0.6, fcb = 20, kc = 0.05, d^^ = 26fi\ and (a) rbr = 3, (b) 
Tbr = 20. The results were obtained according to three different methods, the simplified continuum 
model (bars), the full continuum model (blue solid line) and stochastic simulations (black dashed 
line). 



32 



(ii3j 






- '"'nw 






1? " 


rotated N '\ 

frame ) '\ , 


7 




(a) 




(d) 




(e) 



(f) 



FIG. 4: (a) and (b) Within a rotated coordinate frame, the continuum model (Eq. M) is apphcable 
to treat filament networks growing behind a skewed linear obstacle, (a) Lab frame: A generic skewed 
linear obstacle with skew angle ip. (b) Obstacle frame: The coordinate system is rotated by its skew 
angle if counter-clockwise. In this frame the initial skew of the system is expressed by an additional 
finite lateral motion at velocity v^^ of an otherwise flat obstacle moving at v^^. (c) Parabolic (blue) 
and cosine (red) nonlinear obstacle shapes. The solid lines indicate the nonlinear obstacle shape, 
while the dashed line corresponds to a piecewise-linear triangular approximation. Horizontally 
the parabolic shape has a finite width, while the cosine obstacle shape is analyzed in periodic 
conditions laterally, (d-f ) Piecewise-linear approximations of a parabola-like obstacle geometry, (d) 
The simplest approximation is a flat linear obstacle, (e) Due to the left-right symmetry in obstacle 
shape, in higher order approximations of the nonlinear obstacle (e.g. the triangular approximation 
shown here) , only half of the number of sections have to be considered explicitly. The boundary 
condition at the center of the box has to be adjusted accordingly, (f) Subdividing the piecewise- 
linear sections again and again leads to ever higher accuracy in the approximation of the nonlinear 
obstacle shape and to ever smaller horizontal width of each subsection, d^^ ■ = (i^^/n. Approaching 
the continuum limit, n — )• oo, yields a PDE-mcxiel for nonlinear obstacle shapes. 
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FIG. 5: (a-d) Representative results for the stationary filament orientation distributions of a linear 
tilted obstacle obtained in the PDE- model (thick red solid line), the ODE- model (Eq. (IT])) (blue 
solid line) and in stochastic network simulations (black dashed line). Where applicable also the 
transformation of the resulting PDE-model distributions to the obstacle frame is shown (thin red 
line). Here the orientation patterns can be interpreted as ±35 and +70/0/— 70 patterns, that were 
also active in the tail of flat linear obstacles (cf. Fig. ^ &: Fig. ^. The active parameters are, 
if = 20°, k^ = 20, d^j, = 26m and (a) ry,,. = 3, Unw/'Vfii = 0.8. (b) rbr = 100, Unw/'Vfii = 0.8. 
(c) rbr = 3, fnw/wfii = 0.2. (d) rbr = 100, I'nw/^fii = 0.2. (e) and (f) Comparison of the steady 
state solution for a nonperiodic linear flat obstacle from the ODE-model (Eq. ([l])), the PDE-model 



(Eq. (32)) and stochastic network simulations, (e) For small rbr = 3, the ODE-solution (blue line) 
predicts the absence of a ±35 orientation distribution (cf. Fig. [2]). However, for relatively low 
network velocity at fnw/'Wfii = 0.2, the spatially resolved PDE-solution (red line) as well as network 
simulations (black dashed line) yield such a pattern, (f) The spatial filament distributions are far 
from uniform. Results from PDE-model (solid) and stochastic network simulation (dashed) rather 
show an accumulation of filaments at the lateral flanks of the obstacle, and therefore also horizontal 
outgrowth rates and as a direct consequence, the resulting orientation distribution averaged over 
the whole branching region differs from the prediction of the ODE-model. 
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FIG. 6: Typical stationary filament distributions in the tail of a nonperiodic parabolic obstacle 
(cf. blue solid line in Fig. k[c)) for the parameter combinations, (p = 20°, rbr = 50, d^^ = 2(5fii, 
/cb = 20, kc = 0.05; (1) fnw/i'm = 0.85; (2) Vnw/^m = 0.3. (a) Filament orientation distributions 
spatially averaged laterally over the left hand side half space behind the symmetric obstacle. For 
comparison the results from PDF-model, ODF-model in triangular approximation (as illustrated 
by the dashed blue line in Fig. H[c)) and stochastic simulations are shown, (b) Spatial distributions 
of characteristic filament angles at 6 = 0°, ±35° and ±70°. Interpretation in the lab frame indicates 
coexistence of the competing ±35 and ±70/0/— 70 patterns at different lateral positions along the 
obstacle, (c) Heat map of the spatially resolved filament orientation distributions, where cooler 
colors indicate increasing filament density. Plotting dashed lines at characteristic angles at 9 = 0°, 
±35° and ±70° in the obstacle frame indicates that in this frame a single familiar pattern dominates. 
The apparent coexistence of different patterns is artificially introduced due to a misinterpretation 
of the results in the lab frame. 
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FIG. 7: Same as Fig. |6] for a periodic-cosine shaped obstacle as illustrated by the solid red line 
in Fig. kfc). Active parameters are <f = —20°, rbr = 50, d^^ = 26fi\, k\, = 20, k^ = 0.05; (1) 
f^nw/^fii = 0.5; (2) fnw/^m = 0.1. Again an apparent coexistence of competing patterns disappears 
when results are interpreted in the obstacle frame, that is locally orthogonal to the obstacle surface. 
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